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Atmospheric  Aerosol  Scattering 
Background  Observations 


1.  INTRODUCTION 

Of  particular  importance  to  the  problem  of  electro-optical  (E/O)  propagation  in  the  terrestrial 
atmosphere  is  the  question  of  the  existence  of  atmospheric  aerosols  at  all  altitudes  that  may  affect 
such  propagation.  A  series  of  balloon  flights  with  a  large  polar  nephelometer  provides  results  that 
consistently  indicate  the  presence  of  aerosol  particles,  even  during  periods  relatively  free  of  volcanic 
dust  incursions.  Measurable  scattering  properties  of  aerosols  were  observed  in  both  the  troposphere 
and  stratosphere  and  no  region  was  found  to  be  aerosol-free. 

The  thrust  of  this  report  is  therefore  to  demonstrate  that  ground-based  E/O  systems  operating  in 
the  visible  spectral  range  such  as  LIDARS  will  give  unreliable  results  if  it  is  assumed  that  regions 
identified  by  minima  in  backscatter  contain  no  aerosols.  Similarly,  very  small  particles  (Aitken 
nuclei)  cannot  be  excluded  as  sources  of  error  for  systems  operating  at  wavelengths  below  about 
0.5  pm,  because  they  are  virtually  undetectable  by  LIDAR  alone,  although  they  are  ubiquitous  and  may 
exist  in  large  number  concentrations,  even  in  very  clear  air.  Scattering  from  small  particles  is  not  a 
trivial  problem,  when  molecular  scattering  is  considered  as  a  calibration  threshold  in  remote  aerosol 
or  molecular  density  measurements.  Any  contribution  from  these  particles  results  in  some  error,  if 
the  instrumentation  cannot  distinguish  between  small  particles  and  air  molecules. 

Finally,  two  sets  of  in-situ  data  from  the  balloon-borne  nephelometer  show  contrasting 
scattering  behavior  for  high  concentrations  of  small  particles  in  one  case  and  very  low 
concentrations  of  large  particles  {>  0. 15  pm)  in  the  other,  for  regions  in  the  troposphere  assumed  to  be 
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virtually  aerosol-free  by  some  investigators.1, 2,3,4  Inferences  from  these  results  are  compared  with 
those  obtained  using  particle-counting  methods. 

2.  OBSERVATIONS 

The  balloon-borne  nephelometer  has  been  described  in  previous  publications.5,6’7  The  primary 
components  of  the  Instrument  were  a  150-W  xenon  light  source  and  five  photometers  all  mounted  on 
an  optically  rigid  frame  so  as  to  measure  the  scattered  light  at  15,  30,  50,  100,  and  150°,  from  a  defined 
volume  (about  250  cm1 2 3 4 5 6 7)  of  atmosphere.  Spectral  filters  at  0.475,  0.515,  0.660,  and  0.745  pm  In  the 
photometers  were  used  to  examine  the  properties  of  the  scattered  light  in  the  visible  region.  In 
addition,  the  polarization  at  0.475  and  0.660  pm  was  analyzed  for  the  50,  100,  and  150°  scattering 
angles.  A  series  of  balloon  flights  using  the  equipment  measured  the  altitudinal  variability  in 
scattering  properties  of  atmospheric  aerosols. 

The  optical  parameters  least  affected  by  instrument  calibration  on  the  one  hand  and  most 
sensitive  to  atmospheric  aerosol  scattering  characteristics  on  the  other,  are  the  polarization  of  the 
scattered  light  and  the  dissymmetry  in  angular  scattering  as  a  function  of  altitude.  At  scattering 
angles  near  90°.  the  polarization  is  maximized  by  scattering  from  air  molecules  so  that  a  reduction  in 
polarization  is  attributed  to  atmospheric  aerosols.  Similarly,  the  ratio  of  forward-to-backward 
scattering  or  the  so-called  dissymmetry  parameter,  senses  the  enhanced  forward  scattering  due  to 
aerosols  as  opposed  to  the  characteristic  angular  symmetry  of  molecular  scattering.  Together  with  the 
absolute  scattering  intensity  data,  these  two  parameters  provide  a  means  for  identifying  the  vertical 
structure  of  atmospheric  aerosols  using  optical  diagnostic  techniques. 

Some  experimental  results  from  previous  balloon  flights  have  been  published  as  noted  above. 
The  objective  here  is  to  highlight  certain  aspects  of  these  data,  and  to  include  previously  unpublished 
data,  in  an  attempt  to  provide  evidence  for  the  existence  of  a  background  of  aerosols  in  the 
atmosphere.  These  results  are  pertinent  because  they  were  obtained  during  periods  that  were 
relatively  free  from  volcanic  dust  contamination  of  the  stratosphere  and  might  be  "  msidered 


1.  Northam,  G.B.  et.  al.  ( 1974)  Dustsonde  and  lidar  measurements  of  stratospheric  aerosols:  a 

comparison.  Appl.  Opt.  13:2416-2421. 

2.  Russell,  P.B.  et.  al.  (1975)  Results  of  stratospheric  lidar  observations,  Fourth  Conference  on  the 

Climatic  Impact  Assessment  Program,  Cambridge,  MA.  4-7  February. 

3.  McCormick,  M.P.  and  Fuller,  W.H.  (1975)  Lidar  measurements  of  intense  stratospheric  dust 

layers,  Appl.  Opt.  14:4-5. 

4.  Barteneva,  O.D.  ( 1960)  Scattering  functions  of  light  in  the  atmospheric  boundary  layer,  Bull 

Acad.  Set  USSR  Geophys.  Set,  1237-1244. 

5.  Gibson,  F.W.  and  Dearborn,  F.K.  (1971)  Atmospheric  Optics  Measurements  with  a  Balloon-Borne 

Nephelometer.  AFCRL-TR-7 1-0455,  AD736408. 

6.  Gibson,  F.W.  and  Volz,  F.E.  (1972)  High  altitude  measurements  of  the  optical-scattering 

properties  of  the  atmosphere,  presented  at  Optical  Society  of  America  Meeting,  New  York.  NY, 
April. 

7.  Gibson,  F.W.  (1976)  In-situ  photometric  observations  of  angular  scattering  from  atmospheric 

aerosols,  Appl  Opt  15:2520-2533. 
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representative  of  "normal”  conditions  for  altitudes  up  to  26  km.  Hie  last  balloon  flight  was  made 
from  Holloman  AFB,  New  Mexico  In  June  1976,  approximately  two  years  after  the  incursion  of 
volcanic  dust  from  the  eruption  of  Fuego  in  Guatemala  in  October  1974.8 


2.1  Tuxbidity 

Figure  1  is  a  plot  of  turbidity  as  a  function  of  altitude  and  wavelength.  This  profile  was  obtained 
by  using  the  50°  scattering  intensity  data  to  separate  the  molecular  from  aerosol  scattering.  A  typical 
plot  of  the  intensity,  presented  as  the  volume  scattering  function.  Is  shown  in  Figure  2.  The  latter 
profile,  Pt(50o)  is  due  to  both  aerosol  and  molecular  scattering.  Single  scattering  is  assumed  so  that 
the  intensities  are  separable  and  additive.  The  total  intensity  is 

Pt  =  Pr+  Pa  M 

where  pR  =  volume  scattering  function  due  to  molecular  or  Rayleigh  particles 

{cm-1  sr_1). 

PA  =  volume  scattering  function  due  to  aerosol  particles  (cm-1  sr_1l 
Generally, 

p.,%)  =  oRPR(<t»s)  +  aAPA(^s)  (2) 

where  aR  and  oA  are  respectively  the  Rayleigh  and  aerosol  scattering  coefficients  [cnr1]  and  PR.  PA  are 
the  Rayleigh  and  aerosol  phase  functions  (sr-1).  The  rationale  for  this  formulation  follows  from  the 
empirical  and  theoretical  observations  that  the  Rayleigh  and  aerosol  phase  functions,  at 
approximately  50°  scattering  angle,  are  linearly  related  for  aerosol  size  distributions  found  in  clear 
air.  Barteneva4  has  shown  that  the  phase  function  for  aerosols,  associated  with  atmospheric 
conditions  of  high  visibility,  intersect  In  the  vicinity  of  50°.  Thus 

Pt  _ ,  Pa(5<^  <TaPa(50°)  <tA 

Pr  ~  Pr(50°)  ~  oR  Pa(50°)  -  UoR  ' 

Since  C  =  1 ,  we  obtain  the  turbidity  x,  a  quantity  proportional  to  the  relative  number  of  aerosols  to 
molecules,  given  by 


In  Figure  1  particular  notice  might  be  given  to  the  low  values  of  turbidity  encountered  near  the 
tropopause  at  approximately  10  km  since  this  minimum  is  characteristic  of  LIDAR  observations. 
However  further  comments  will  be  made  after  a  more  detailed  analysis  is  presented  later. 


8.  Meinel,  A.B.  and  Meinel,  M.P.  (1975)  A  stratospheric  dust  event  of  November  1974,  Science 
188:477-481. 
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An  earlier  flight  In  November,  19705  resulted  In  the  turbidity  profile  shown  in  Figure  3. 
Although  the  values  found  in  the  stratosphere  (designated  the  Junge  aerosol  layer®)  were  very  close  to 
those  observed  in  1976,  the  tropospheric  values  were  strikingly  different.  An  explanation  for  the 
differences  is  given  later  in  this  report  and  only  the  close  comparison  of  the  stratospheric  data, 
particularly  at  0.660  pm  wavelength,  is  emphasized  at  this  Juncture.  The  latter  comparison  is 
particularly  meaningful  because  of  the  duration  of  about  six  years  between  flights;  a  period  that 
Included  a  major  volcanic  eruption;  that  is,  the  Fuego  event.8 


2.2  Polarization 


Figure  4  is  a  plot  of  the  polarization  ratio  for  the  1976  flight.  This  profile  was  determined  by 
taking  the  ratio  of  100°  scattered  light  at  0.475  pm.  measured  sequentially  through  the  filters 
containing  linear  polarizers.  The  polarizers  were  oriented  so  as  to  detect  the  vertical  and  horizontal 
components  of  the  scattered  light.  Figure  5  is  a  similar  plot  for  the  1970  data,  but  the  photometer  was 
located  at  the  90°  scattering  angle.  In  these  profiles,  the  parameter  displayed  is  related  to  the  degree  of 
linear  polarization  by 


Pol. 


_  ly~  4i  _  WIh  1 
lv  +  Ih  fy/f  h+1 


(5) 


where  Iv  and  IH  are  the  scattered  Intensities  perpendicular  and  parallel  to  the  plane  of  scattering  and 
the  quantity  Iv/In  is  an  absolute  measurement,  because,  for  a  specific  photometer/ detector,  the 
calibration  factor  is  common  to  Iv  and  I|{.  Further  analysis  of  this  parameter  has  been  given  In 
previous  publications7  and  the  profiles  are  presented  here  primarily  to  illustrate  differences  in  the 
observations  for  the  two  flights  as  well  as  to  indicate  the  consistency  between  the  polarization  ratio 
and  other  data,  such  as  turbidity  profiles.  In  terms  of  comparison,  the  1976  data  indicate  the  presence 
of  large,  so  called,  Mie9  particles,  at  virtually  all  altitudes  inasmuch  as  the  polarization  ratio  reaches 
peak  values  between  7  and  8  in  the  9  -  10  km  region,  which  coincides  with  the  turbidity  minimum. 

This  maximum  corresponds  to  76  percent  polarization  of  the  scattered  light.  An  aerosol-free  region 
requires  a  ratio  of  about  14  (87  percent  polarization)  which  would  indicate  pure  molecular  scattering. 
The  1970  data  in  Figure  5  show  maxima  between  9  and  10  at  several  points  in  the  5-10  km  altitude 
range,  indicating  about  82  percent  polarization,  while  molecular  scattering  would  give  an  Iv/Iu  value 
of  24.  that  is.  92  percent  polarization  for  the  90°  scattering  angle.  In  both  profiles  however,  the  ratios 
are  approximately  the  same  for  the  Junge  layer,  above  16  km,  when  the  10°  difference  in  photometer 
location  is  considered.  That  is,  a  comparable  reduction  in  polarization  occurred  when  the 
stratospheric  Mie  particles  were  encountered.  Noteworthy  also  is  that  such  reduction  (to  virtually  the 
same  value)  also  occurred  in  the  troposphere  for  the  1976  data,  thereby  suggesting  large  particles  were 
present,  which  affected  the  polarization  in  much  the  same  manner  as  those  in  the  stratosphere. 


9.  Junge.  C.E.,  Chagnon,  C.W.,  and  Manson,  J.E.  (1961)  Stratospheric  aerosols,  J.  Meteorol  18: 
81-98. 
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Figure  3.  Turbidity  cta/Or  vs  Altitude  (November  1970) 


4  5  6  7  8  9  10  II 

4750A(V)/4750A(H)  (90°  PHOTOMETER) 

Figure  5.  90°  Polarization  Ratio  vs  Altitude  (November  1970) 
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2.3  Dissymmetry 

It  has  been  shown7  that  the  ratio  of  scattering  intensities  for  two  symmetric  scattering  angles 
can  be  related  to  the  turbidity  as  a  function  of  altitude  and  as  a  result  a  parameter  is  obtained  that 
indicates  the  altitudinal  variability  in  angular  scattering  by  aerosol  particles.  The  symmetric  angles 
utilized  in  the  nephelometer  flights  were  30°  and  1 50°  and  thus  the  dissymmetry  is  defined  as: 

Px(30°)  PKX  (30°)  +  (oA/oR)  PAX  (30°) 

R  “  ^(150°)  "  PRX  (150°)  +  (oA/oR)  PAX  (150°)  1  ' 

where 

R  s  dissymmetry  factor  or  simply  the  forward-to-backscatter  ratio, 

PRX  =  normalized  Rayleigh  phase  function, 

PAX  =  normalized  aerosol  phase  function,  and 

oA/oR  =  x.  the  turbidity  (a  quantity  independent  of  scattering  angle). 

Now  PrA.  (30°)  =  PRX  (150°)  because  Rayleigh  scattering  is  symmetric  relative  to  the  90°  scattering  angle. 
Aerosol  particles  with  a  size  parameter,  X  =  0.8,  scatter  more  light  in  the  forward  direction10 

which  means  PAX  (30°)  >  PaA.  (150°).  Consequently.  R  >  1  when  x  >  0,  with  the  result  that  the  presence  of 
aerosols  should  be  evidenced  by  values  of  the  R-parameter  greater  than  unity.  Moreover,  this  factor  is 
wavelength  dependent,  because  the  turbidity  varies  with  wavelength,  and  it  will  be  shown  in  the  next 
section  that  such  wavelength  dispersion  provides  quantitative  information  on  the  aerosol  phase 
function. 

The  vertical  profile  of  R  for  the  1976  flight  is  shown  in  Figure  6.  This  plot  is  simply  the  ratio  of 
the  experimentally  determined  scattering  functions,  p\  (30“)/pX  (150°),  at  the  30  and  150“  angles,  and 
the  0.475.  0.515,  and  0.660  pm  values  are  displayed.  This  profile  is  clearly  very  similar  in  overall 
structure  to  the  turbidity  plot  shown  in  Figure  1  suggesting  a  correlation  between  R  and  x  as  predicted 
by  Eq.  (6).  The  Junge  layer  is  delineated  in  particular  by  the  R-parameter  with  peak  values  occurring 
at  about  18  km,  and  prominent  dispersion  toward  the  red,  (0.660  pm). 


3.  DISCUSSION  OF  THE  OBSERVATIONS 
3.1  Theoretical  Considerations 

More  quantitative  interpretations  of  the  data  from  the  balloon-borne  nephelometer 
experiments  require  an  outline  of  applicable  theory.  These  concepts  are  either  extensions  of  or  direct 
use  of  Rayleigh  and  Mie  scattering  theory,  which  have  served  to  explain  atmospheric  light  scattering 
phenomena  generally. 


10.  Van  de  Hulst,  H.C.  (1957)  Light  Scattering  by  Small  Particles,  John  Wiley  &  Sons,  New  York. 
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In  view  of  the  fact  that  very  small  atmospheric  particles  cover  a  size  range  from  air  molecules 
Ur  3  x  1CT4  pm)  to  condensation  nuclei  or  Aitken  particles  (=10-2  pm)  It  Is  necessary  to  examine 
Rayleigh  scattering  as  a  function  of  particle  size.  Van  de  Hulst10  has  shown  that  when  the  particle  size 
parameter. 


X  = 


2nr 


<  0.8, 


the  angular  scattering  pattern  Is  virtually  symmetrical.  Recent  laboratory  studies  with  laser  light 
have  also  verified  that  negligible  error  will  be  Introduced  when  this  upper  limit  Is  assumed  for  the  size 
parameter  of  Rayleigh  scatterers. 1 1  The  Implication  for  the  balloon-borne  nephelometer  Is  that  at 
X  =  0.475  pm,  atmospheric  particles  with  radii  less  than  0.06  pm  can  be  considered  to  be  Rayleigh 
scatterers.  This  limit  therefore  Includes  Aitken  nuclei,  which  are  apparently  critical  to  cloud 
formation  and  have  been  observed  In  the  troposphere  with  number  concentrations  the  order  of  104 
cm-3.  The  Rayleigh  scattering  formula  Is  usually12  expressed  as 

fR&t>.  X)  =  7£it2  tVr^2  ~  2  J  (1  +  ros2  I*311"1  -  sr_1l  (7) 

where 

fR  =  angular  volume  scattering  function 
nR  =  number  particles  cm3 
m  =  refractive  index 
v  =  volume  of  particles 
X  =  wavelength  of  scattered  light 
<t>  =  scattering  angle 

For  spherical  particles  (such  as,  air  molecules)  this  equation  becomes 

fR(<fr,  X)  =  8tc4^”2  +  ^J  n*jj-  (1  +  cos2  <)>)  [cm-1  -  sr"1]  (8) 

where  r  =  radius  of  particles,  assumed  to  be  monodlspersed.  The  scattering  function  is  therefore 
proportional  to  the  6th  power  of  particle  radius.  Indeed,  Lord  Rayleigh  originally  determined  the  X-4 
wavelength  dependence  from  dimensional  considerations.  That  Is,  the  scattering  Intensity  is 
proportional  to  the  square  of  the  particle  volume  and  the  number  density,  and  the  X-dependence  of  X-4 
Is  required  for  the  Intensity  to  have  the  proper  dimensions. 


11.  Sassen,  K.  (1981)  Infrared  (10.6  pm)  scattering  and  extinction  In  laboratory  water  and  Ice  clouds, 

Appl  Opt  20: 185-193. 

12.  Kerker,  M.  (1969)  The  Scattering  of  Light  and  Other  Electromagnetic  Radiation,  Academic  Press, 

New  York. 
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Because  one  of  the  objectives  In  this  report  is  to  emphasize  the  Importance  of  the  r6  factor  In 
Rayleigh  scattering  formulations,  It  is  Informative  to  illustrate  that  this  relationship  can  also  be 
deduced  from  Mle  scattering  considerations.  Bullrich13  has  examined  the  scattering  properties  of 
atmospheric  aerosols  that  have  number  densities  which  follow  a  power-law  type  distribution  in  size. 
Such  distributions  were  found  by  Junge9  in  balloon-borne  investigation5  of  naturally  occurring 
aerosols  and  are  given  by 


dn 

d(log  r) 


Cr_v  Icrrr3] 


(9) 


This  distribution  gives  the  number  of  particles  dn  with  radius  between  r  and  (r  +  dr).  This  can  be 
rewritten  in  the  form 

^  =  0.434  Cr-<v+1>  Icrrr4]  . 

where  the  change  of  variable  is.  d(Iog  r)  =  -p-  =  0.434  — . 

For  such  a  size  distribution  Bullrich  determines  the  scattering  function  to  be 

*2 


(10) 


fAW>.  =  0.434 


§(¥T  J 


lf».  X) 

yVt  1 


dx  [cm-1  -  sr"1] 


(11) 


where  x  =  ,  is  the  size  parameter.  The  integral 


4 

/ 


1(4>.  X) 

TrV+l 


dx  =  T](4>,  X) 


(12) 


Is  the  angle  dependent  factor  or  phase  function  that  gives  the  angular  scattering  behavior  in  terms  of 
the  size  parameter,  which  is  dimensionless.  This  factor  T]fo>,  X)  can  be  calculated  for  a  given  size 
distribution  and  wavelength,  using  Mie  scattering  computer  codes.  The  constant,  C,  is  determined 
from  the  total  number  of  particles  In  a  particular  size  range.  That  is,  the  total  number  N.  of  particles 
between  rj  and  r2  is 

? 

N=  J  dn  =  0.434 C  J  Hv+1>  dr,  from  Eq.  (10). 


Integration  gives 


0.434  C 
vr,v 


(13) 


13.  Bullrich,  K.  (1964)  Scattered  radiation  in  the  atmosphere  and  the  natural  aerosol,  Advances  in 
Geophysics  10:101-257. 
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If  rj  «  r2,  then  N =0.434  C 


- and  for  some  r,  and  rA, 

v  r,v  l  a 


NvrAv 

C  ~  0.434  ' 


(14) 


When  this  value  is  used  in  Eq.  (1 1)  we  obtain 
NT  ( \r^ 

fA(0.x)  -yItJ  vrAVl)^'  -  {15) 

For  v  =  6  we  get 

fA(0.  X)  =  3N  (^)  4  rA6  T1(0,  X)  =  3N  ^  rA6  rfo.  X)  . 

The  explicit  dependence  of  fA,  on  X-4  and  rA6  is  clear.  Although  Eq.  (15)  will  be  applied  later  in  the 
analysis  of  the  turbidity  data,  the  special  case  v  =  6  simply  illustrates  the  consistency  between 
Rayleigh  and  Mie  theoiy  in  the  domain  of  small  particles.  It  should  be  noted  however  that  the  radius 
exponent  and  the  wavelength  exponent  are  both  functions  of  v.  That  is.  from  Eq.  (15). 

fA  ~  X-v  +  2  =  XA 

where  a  =  v  -  2 

and  v  is  the  exponent  in  the  power-law  size  distribution.  This  wavelength  factor  is  the  well-known 
Angstrom  relation  and  has  undergone  extensive  empirical  verification,  particularly  under  conditions 
where  an  exponential  size  distribution  was  apparent. 


3.2  A  Postulated  Size  Distribution 

Because  optical  scattering  effects  of  very  small  aerosols  are  usually  neglected,  since  they  have 
little  effect  on  visibility,  few  experimental  data  exist  relevant  to  their  influence  on  atmospheric 
scattering  when  large  particles  are  virtually  absent,  that  is,  under  very  clear  sky  conditions.  To  assess 
such  an  environment  we  constructed  a  particle  size  distribution  with  a  dearth  of  particles  greater  than 
0. 1  pm  in  radius. 

An  inverse  power-law  distribution  of  the  type  suggested  in  Eq.  (9)  with  the  exponent  v  =  6.  implies 
a  constant  slope  or  drop-off  in  particle  radius  from  about  0.04  pm  to  approximately  0.07  pm.  However, 
actual  measurements  of  Aitken  particles  by  Junge14  and  others15,16,  have  indicated  that  such 
particles  exhibit  a  log  normal  size  distribution,  probably  due  to  coagulation  of  particles  with  radius 
less  than  about  0.02  pm.  Moreover,  particles  with  radii  greater  than  this  lower  limit  comprise  the 


14  Junge,  C.E.  (1961)  Vertical  profiles  of  condensation  nuclei  in  the  stratosphere,  J.  Meteorol. 
18:501-509. 

15.  Went,  F.W.  (1964)  The  nature  of  Aitkin  condensation  nuclei  in  the  atmosphere.  Nature  51: 1259- 

1267. 

16.  Whitby,  K.T.  (1975)  Modeling  of  Atmospheric  Aerosol  Particle  Size  Distribution,  Prog.  Rep.  253, 

Particle  Technology  Laboratory.  Univ.  of  Minnesota,  Minneapolis. 
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bulk  of  the  active  condensation  nuclei  in  the  formation  of  clouds.17  Investigations17  into  the  process 
of  water  vapor  condensation  involve  the  so-called  Kelvin  factor,  which  indicates  that  condensation 
(that  is,  clouds,  fog,  etc)  cannot  begin  on  particles  less  than  about  0.02  pm  radius  due  to  an  abnormally 
high  degree  of  supersaturation  required.  Under  normal  meteorological  conditions  in  the  atmosphere, 
hygroscopic  nuclei  serve  as  the  centers  on  which  water  vapor  condenses.  Indeed,  several  studies  have 
established  correlations  between  spectra  of  cloud  particles  and  those  of  small  aerosol  particles  from 
which  they  evolve,  depending  upon  factors  such  as  relative  humidity  and  supersaturation.  In  this 
connection  Nilsson18  has  used  a  log  normal  size  distribution  given  by 


dn  2.3Nt  f  1  flnr/rm^21 

"  V2rc"  In  SG  i13  ”  2  V  In  SG  )  \ 


(16) 


where  ^  ~py  =  number  density  per  log  radius  interval  [cm_3l 
Nt  =  total  number  density  [cm-3] 
rm  =  median  particle  radius  (pm) 

In  SG  =  standard  deviation  of  In  r. 


To  determine  an  appropriate  distribution  for  the  size  range  of  small  particles  being  considered  it  is 
only  necessary  to  define  In  SG  so  that  most  of  the  particles  will  have  radii  in  the  range  0.02  pm  <  r  < 
0.07  pm.  This  can  be  accomplished  by  matching  the  power-law  distribution  which  dropped  ofT  as  r-6, 
and  the  log  normal  distribution  at  an  initial  value  and  at  r  =  0.07  pm.  That  is,  the  respective 
distributions  have  the  form 


and 


dn  _ 
d(log  r)  -  nro 

dn  _ 
d(log  r)  ~  nr° 


(r)6 

[“P-HtStKt)2]-  wheren 


2.3Nt 
A  2  In  SG 


If  we  then  set 


1  ( Inr/ro 

=  “P  "  2  VETST 


) 


2 


we  obtain 


6  In  (r/ro)  = 


1  flnr/rol2 

2  l  In  SG  ) 


and 


InSG  =  0.289  Vln (r/ro)  . 
I 


(17) 

(18) 

(19) 

(20) 

(21) 


17.  Junge,  C.F.  and  McLaren,  E.  (1971)  Relationship  of  cloud  nuclei  spectra  to  aerosol  size 

distribution  and  composition,  J.  Atmos.  Set  28:382-390. 

18.  Nilsson,  B.  (1979)  Meteorological  influence  on  aerosol  extinction  in  the  0.2  pm  -  40  pm 

wavelength  range,  Appl.  Opt  18:3457-3473. 
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Now  assuming  r  <  0.07  pm  and  considering  r0  =  0.04  pm  as  the  median  radius  we  find 
In  SG  =  0.216 

as  the  standard  deviation,  or  the  spread  of  our  distribution,  with  a  2(fn  SG)  value  (containing  95 
percent  of  the  particles)  of 

0.026  pm  <  r  <  0.062  pm . 

A  representative  value  or  nro  from  the  measurements  of  Junge  is 
n^,  =  5  x  lOVcm3  at  rD  =  0.04  pm . 

The  tropospheric  model  suggested  by  Shettle  and  Fenn19  [converted  to  dn/d(log  r)|  has  nro  =  1.07  x 
104/cm3  at  r  =  0.04  pm  for  a  log  normal  distribution  with  o  =  0.35  and  NT  =  lOVcm3.  The  distribution 
suggested  here  gives  Nj  =  1. 17  x  104/cm3,  for  =  5  x  104/cm3.  The  total  number  of  particles  Is  thus 
close  to  those  In  other  models,  but  the  size  distribution  Is  restricted  to  particles  less  than  0. 1  pm  In 
radius.  Such  a  restriction  does  not  seem  unreasonable  for  extremely  clear  air,  such  as  that  found  In 
the  desert  environment  of  New  Mexico. 

The  suggested  size  distribution  was  used  to  determine  the  pertinent  scattering  parameters  using 
AFGL  Mie  codes.  The  angular  symmetry  and  high  degree  of  polarization  at  90°  demonstrated  a 
scattering  pattern  with  distinctive  Rayleigh  characteristics.  The  computed  values  of  volume 
scattering  functions  were: 

1.94  x  KT4  km-1  sr-1  at  0.475  pm 
1.42  x  10-4  km-1  sr-1  at  0.515  pm 
5.30  x  1CF5  km-1  sr-1  at  0.660  pm. 

These  results  clearly  Indicate  a  X-4  dependence  upon  wavelength  within  less  than  2  percent  error. 
Finally,  this  distribution  of  small  particles  would  give  an  expected  value  of  turbidity  of  0.34  at  the  10 
km  altitude  level  often  cited  as  a  region  of  minimum  scattering.  However,  these  particles  alone  could 
not  account  for  the  abnormal  dispersion  in  turbidity  observed  In  1970  because  of  the  X-4  dependence. 
Nonetheless,  the  resulting  contribution  to  the  total  scattering  cannot  be  assumed  to  be  Insignificant 
since  It  would  introduce  an  error  of  at  least  30  percent  if  this  region  were  considered  aerosol  free.  For 
the  nephelometer  data,  such  high  cencentrations  of  small  particles  provide  an  additive  factor  that  Is 
examined  qualititively  in  the  next  section  to  show  that  under  these  conditions  the  normal  Angstrom 
relation  Is  modified. 


3.3  Small  Particle  Effects  on  the  Observations 

The  postulated  distribution  has  virtually  no  particles  with  radii  greater  than  0. 1  pm.  A  more 
realistic  picture  of  atmospheric  aerosol  particles  would  Include  some  particles  >0.1  pm  In  radius. 


19.  Shettle,  E.P.  and  Fenn,  R.  W.  (1979)  Models  for  the  Aerosol  of  the  Lower  Atmosphere  and  the 
Effects  of  Humidity  Variations  on  their  Optical  Properties,  AFGL-TR-79-0214,  ADA085951. 
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particularly  In  view  of  the  well-documented  enhanced  concentration  of  large  particles  above  the 
tropopause.  that  comprises  the  Junge  layer  In  the  stratosphere.  It  is  therefore  reasonable  to  expect  a 
small  number  of  these  aerosols  to  be  found  below  the  tropopause  due  to  atmospheric  dynamics  as  well 
as  settling  processes.  Because  the  sources  of  small  particles  and  large  particles  are  assumed 
uncorrelated  and  independent,  the  size  distribution  suggested  here  is  bimodal.  Other 
investigators20-21  were  compelled  by  their  experimental  data  on  sky  radiation  and  absorption  to 
postulate  the  existence  of  steep  size  distributions  of  small  particles  along  with  independent 
distributions  of  large  particles.  That  a  comparable  situation  existed  in  the  troposphere  during  the 
1970  balloon-borne  nephelometer  flight  was  suggested  as  an  explanation  for  the  "unusual"  dispersion 
that  resulted  when  the  turbidity  was  plotted  as  a  function  of  wavelength  as  seen  in  Figure  3. 

The  basic  approach  in  analyzing  the  1970  turbidity  data  was  to  assume  that  the  turbidity 
emanated  from  non-molecular  Rayleigh  particles  and  Mie  particles,  that  is,  those  with  radii  >0.1  pm. 
Because  the  scattering  intensities  are  separable,  the  observed  volume  scattering  function  for  a 
bimodal  distribution  of  aerosol  is, 

PAobs=Pp  +  PA  (22) 

where 

Pp  =  volume  scattering  function  for  Rayleigh  particles 

PA  =  volume  scattering  function  for  Mie  particles 
the  observed  turbidity  is  then, 

=  1231 

where  PR  =  volume  scattering  function  for  air  molecules. 

When  this  expression  is  rewritten  as  a  product  of  two  factors 

1241 

we  see  that  when  PP  >  0,  the  measured  turbidity  is  greater  than  that  due  to  Mie  particles  alone  by  the 
factor  in  brackets.  Moreover,  the  ratio  (Pp/PA)  provides  an  additive  term  that  indicates  the  relative 
contribution  of  small  particles  versus  large  particles.  In  principle  then,  if  pp  is  greater  than  PA  this 
term  is  indeed  not  negligible.  To  obtain  a  qualitative  estimate  of  the  magnitude  of  this  factor  consider 
the  following: 

From  Eq.  (15)  an  exponential  distribution  of  Mie  particles  with  v  =  4  would  give  a  scattering 
function  proportional  to  X-2.  A  similar  power-law  distribution  of  Rayleigh  particles  has  v=  6  and  a 
X'4  wavelength  factor.  (Note  that  this  distribution  is  the  same  as  the  log-normal  distribution  in 


20.  Rossler,  F.  ( 1972)  Aerosol  layers  in  the  atmosphere,  Space  Research  XII  29:423-43 1 . 

21.  Pilipowsky,  S.  et.  al.  (1968)  Investigation  of  the  stratospheric  aerosol  by  infrared  and  Lidar 

techniques,  J.  Geophys.  Res.  73:7553-7560. 
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Section  3.2  for  small  particles).  If  the  index  of  refraction  is  approximately  the  same  for  both 
distributions  [and  at  =  50°  the  phase  functions  are  nearly  equal,  see  Eq.  (3)1  we  have 

PL_(2jd2  6  Np  44-59-^f^ 
pA  "  ^  4NAf^F  mN  1  - 


(25) 


when 


NP  =  number  or  Rayleigh  particles  with  radius  rP 


Na  =  number  of  Mie  particles  with  radius  rA. 

Now  if  rp  =  0.05  pm  and  rA  =  0.2  pm  and  X  =  0.475  pm  we  obtain 
Bp  NP 

21=^(2.56x10-3). 


(26) 


Pa 

For  jj£->  2,  it  is  only  necessaiy  that  jr->  103,  with  the  implication  that  103  Rayleigh  particles  are 
Pa  Pa 

required  to  give  twice  the  scattering  intensity  of  1  Mie  particle.  The  next  question  is  how  does  such  a 
factor  (Pp/Pa)  >  2  affect  the  turbidity-wavelength  dependence?  Since  the  balloon-borne  nephelometer 
provided  simultaneous  measurements  at  different  wavelengths,  Eq.  (24)  can  be  employed  to  determine 
the  turbidity  quotient  for  two  wavelengths.  In  particular  it  is  desirable  to  establish  a  condition  under 
which  the  observed  turbidity  at  0.475  pm  is  greater  than  that  observed  at  0.660  pm.  A  ratio  of  such 
turbidities  from  Eq.  (24)  is 


(Pa/Pr)  obs  ■  X,  _  f(PA/pR)  Mie  Xfl  f  1  +  (pP/pA)  Xf 
(Pa/Pr)  °Ps  '  ^2  _ (Pa/Pr)  Mie  X^  _1  +  (Pp/pA)  Xg_ 


(27) 


Now  we  have  seen  previously  that  PA  =  (constant)  X"2  when  v  =  4,  where  the  constant  factor  is 
independent  of  X  even  though  it  is  a  function  of  the  particle  size  distribution.  Moreover,  at  a 
designated  altitude  and  molecular  number  density.  PR  =  (constant)  X-4. 

Thus  (Pa/Pr)  Mie  =  (CA/CR)  X2  (28) 

provides  the  X-dependence  for  the  factor  in  the  turbidity  due  to  Mie  particles.  Since  the  same  particles 
are  scattering  light  at  both  wavelengths  we  get 

(Pa/Pr)  Mle  _  AiV  f9q> 

(PA/pR)MieX2  "(XaJ 

and  Eq.  (27)  becomes 


(Pa/Pr)  Obs  _  Ai'ff  1  +  (Pp/PA)  Xt 
(Pa/PrI  ^2  J  1  +  (Pp/Pa)  ^2 


This  quotient  is  greater  than  unity  when 


1  +  (Pp/P^  X; 
1  +  (Pp/Pa)  ^2 


=  1.93,  for  X!  =  0.475  pmm  and  =  0.660  pm . 


(31) 
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Further  reduction  of  this  expression  leads  to 
(pp/pA)  A.j  >  0.93  +  1.93  (Pp/PaJ  \2  ■ 


(32) 


This  implies  (0p/pA)  0.475  4  >  1.93,  when  (Pp/Pa)  0.660  4  <  =  0.518 . 

This  condition  is  virtually  the  same  as  that  suggested  above  when  the  number  concentration  of  small 
particles  is  greater  than  that  of  large  particles  by  a  factor  of  103.  While  in  principle  this  example 
simply  illustrates  that  the  factor  ( 1  +  Pp/PAI  in  Eq.  (24)  due  to  Rayleigh  particles  has  a  wavelength 
dependence  of  approximately  A.-2,  the  number  of  small  particles  relative  to  large  particles  is  the 
dominant  physical  quantity  affecting  (Pp/PA).  Consequently  the  instrumental  capability  to  detect 
small  particles  is  paramount.  This  situation  can  be  surmised  from  Eq.  (30)  when  (Pp/Pa)  0.660  4  «  1 . 
which  is  perhaps  realistic  for  the  1970  nephelometer  data  as  well  as  lidar  systems  operating  at  red 
wavelengths.  In  that  case  the  fight  hand  side  of  Eq.  (30)  is  about  2  and  is  indeed  consistent  with  the 
turbidity-wavelength  profile  in  the  troposphere  in  Figure  3. 

Even  though  the  size  distribution  of  Mle  particles  was  chosen  for  demonstration  (and  is  not 
unique)  the  results  indicate  that  high  concentrations  of  small  particles  together  with  large  particles 
cause  a  predictable  effect  on  the  turbidity-wavelength  relationship.  Therefore,  the  foregoing  analysis 
appears  to  be  a  plausible  explanation  for  the  "unusual"  dispersion  observed.  Furthermore,  such  an 
explanation  is  consistent  with  the  high  degree  of  polarization  of  the  scattered  light,  also  detected  in 
the  same  altitude  range. 

3.4  Inferences  from  Angular  Scattering  Behavior 

The  observations  from  the  flight  in  1976  were  examined  in  an  effort  to  correlate  the  turbidity 
with  altitudinal  variations  in  angular  scattering.  The  1970  vertical  profiles  of  forward-to-backward 
scattering  were  published  previously5  and  showed  similarities  to  the  1976  data,  despite  the  observed 
turbidity  dispersion.  The  method  of  analysis  presented  here  is  therefore  applicable  to  the  earlier  data, 
but  is  complicated  by  the  problem  of  separating  the  contribution  to  the  turbidity  from  Rayleigh 
particles,  because  in  principle  only  Mle  particles  cause  changes  in  angular  scattering  symmetry.  It  is 
shown,  however,  in  the  next  section  that  the  similarities  between  these  sets  of  data  provide  additional 
evidence  for  the  reliability  of  the  1970  turbidity  profiles,  but  the  1976  data  offer  a  less  complicated 
basis  for  considering  Mie  particle  effects.  Eq.  (6)  contains  the  Rayleigh  phase  function  PR.  the  aerosol 
phase  function,  PA  and  turbidity  x.  The  1976  profile  of  x  did  not  indicate  "unusual"  behavior  in  the 
troposphere  like  that  experienced  in  1970.  That  Is,  the  wavelength  dispersion  is  toward  the  red 
wavelength  as  predicted  by  Eq.  (30)  when  (Pp/PA)  «  1  for  all  wavelengths.  PR  in  Eq.  (6)  therefore  refers 
to  molecular  scattering  and  is  usually  expressed  as 

PR 

PR(<S  =  ~  (1  +  cos2«t>)  IsrM,  and  j  (4>)dn  =  1  (34) 

1671  4 1 
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where  dft  =  the  element  of  solid  angle.  To  facilitate  comparison  with  published  aerosol  phase 
functions  the  normalized  form  used  by  Elterman22  is  more  convenient,  because  It  approaches  unity  at 
=  50°.  That  is.  PR  =  3/4  (1  +  cos2  <fr)  and  PR  (30°)  =  PR(150°)  =  1.3.  Eq.  (6)  then  becomes. 


1.3  +  x  PA(30°) 
R_  1.3  +  x  PA(150°) 


(35) 


where  x  =  — “  is  the  turbidity  determined  from  50°  Intensity  observations.  Rearranging  this  expression 

Or  • 

gives 


Pa(30°)  =  1.3  — — —  +RPa(150°)  . 


(36) 


For  a  given  size  distribution  of  aerosol  particles,  the  phase  function  changes  very  little  with 
wavelength  In  the  visible  spectral  range  and  to  a  first  approximation  PA(X))  =  PA(X^).  For  example  the 
scattering  function  for  the  normal  power-law  particle  size  distribution  was  shown  previously  (see  Eq. 
(15))  to  contain  a  factor  with  the  explicit  wavelength  dependence  and  the  phase  function,  which  is 
dimensionless.  The  phase  function  for  the  two  wavelengths,  Xj  and  X^,  at  a  fixed  angle,  can  differ  by  at 
most  a  small  constant  factor.  While  this  factor  may  be  significant  In  theory,  it  is  virtually 
undetectable  experimentally  in  view  of  the  variability  in  the  observed  quantities  R  and  x  as  a  function 
of  altitude  and  wavelength.  With  this  assumption  that  PA(Xj)  =  P^A^),  simultaneous  measurements  of 
R  and  x  at  X,  and  A^  provide  two  equations  for  PA(150°)  in  the  quantity  PA(150°)  which  can  be  solved  to 
give 

1371 

where  Ri  =  value  of  R  at  Xp  Rj  =  value  of  R  at  X2  and  similar  designations  for  the  measured  values  of  x. 
This  result,  when  used  in  Eq.  (36)  with  either  set  of  values  for  R  and  x,  indicates  that  PA(30°)  and 
PA(150°)  are  determinable  from  the  observed  parameters.  Interestingly.  Eq.  (37)  places  bounds  on  the 
experimental  data  because,  for  PA(150°)  >  0, 

R]  <  ^(R2 -  1)  +  1,  when  R,  > Rj  and  x,  >  .  (38) 


This  method  was  used  to  compute  PA(30°)  and  PA(  1 50°)  at  several  altitudes  from  the  vertical  profiles  of 
x  and  R,  to  arrive  at  experimentally  inferred  angular  scattering  characteristics  of  aerosol  particles 
encountered  during  the  1976  flight.  The  values  are  listed  in  Table  1.  Because  only  two  symmetric 
scattering  angles  were  used,  it  is  ambiguous  and  speculative  to  interpret  the  results  in  terms  of  well 
defined  particle  size  distributions.  Nevertheless,  it  is  possible  to  make  comparisons  with  theoretical 
models  used  by  others. 


22.  Elterman,  L.  (1966)  Aerosol  measurements  in  the  troposphere  and  stratosphere,  AppL  Opt 
5:1769-1775. 
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Table  1.  Scattering  Properties  and  Particle  Number  Concentrations  (June  1976) 


Altitude 

(km) 

Rl 

*1 

R2 

PA(30°) 

PA(  1 50°) 

Pa(30°)/Pa(1  50°) 

V*2 

D 

N(>0. 15mm) 

— 

cm3 

8 

2.2 

0.69 

1.5 

0.24 

3.60 

0.634 

5.73 

2.73 

2.8 

1.23 

2.35 

9 

1.8 

0.30 

1.3 

0.10 

5.03 

0.867 

5.80 

3.0 

2.6 

0.36 

0.77 

10 

2.0 

0.40 

1.4 

0.14 

4.88 

0.833 

5.85 

2.86 

2.8 

0.52 

1.24 

11 

3.3 

3.6 

2.2 

1.2 

2.22 

0.422 

5.26 

3.0 

2.6 

3.4 

9.09 

12 

2.2 

0.9 

1.3 

ESSffi 

2.26 

0.24 

9.42 

4.0 

2.0 

0.61 

1.84 

13 

2.0 

1.0 

1.3 

2.04 

0.37 

5.51 

4.0 

2.0 

0.51 

1.97 

14 

2.1 

1.11 

1.5 

2.92 

0.78 

3.74 

3.0 

2.6 

0.66 

2.55 

15 

2.5 

1.47 

1.5 

0.43 

1.79 

0.185 

9.67 

3.41 

2.3 

0.56 

2.46 

16 

3.2 

2.1 

1.7 

0.50 

2.34 

0.305 

7.67 

4.02 

2.0 

0.76 

3.94 

17 

4.0 

3.8 

2.1 

1.0 

1.87 

0.212 

8.84 

3.8 

2.0 

1.18 

7.30 

18 

4.18 

3.5 

2.2 

1.01 

1.94 

0.182 

10.67 

3.46 

2.2 

0.83 

6.11 

19 

4.7 

3.7 

2.3 

1.10 

1.76 

0.098 

17.96 

3.36 

2.3 

0.75 

6.57 

20 

4.4 

3.2 

2.6 

1.20 

2.24 

0.196 

11.43 

2.67 

3.0 

0.97 

10.08 

21 

3.2 

3.8 

2.0 

1.30 

1.41 

0.206 

6.84 

2.92 

2.7 

0.98 

12.07 

22 

3.05 

3.17 

2.2 

1.50 

1.55 

0.234 

6.62 

2.11 

3.7 

0.88 

12.81 

23 

4.4 

5.0 

3.0 

1.60 

3.21 

0.529 

6.07 

3.13 

2.5 

0.70 

12.03 

24 

2.5 

2.8 

2.1 

1.80 

1.31 

0.245 

5.35 

1.56 

4.6 

0.56 

11.34 

25 

3.5 

5.3 

2.5 

2.0 

1.88 

0.362 

5.19 

2.65 

0.70 

16.70 

26 

2.6 

4.9 

1.8 

2.0 

0.729 

0.119 

6.17 

2.5 

0.60 

16.80 

Since  the  explicit  relationship  between  scattering  intensity  and  wavelength  for  a  power-law 
particle  size  distribution  provides  a  means  to  estimate  particle  number  density.  Eq.  (15)  was  again 
employed  to  examine  the  turbidity-wavelength  data.  By  rearrangement  and  change  of  notation,  Eq. 
(15)  takes  the  form 

n  -  v  (t)2?*!^  |mr31  ■  °9! 

The  factor  q(4>,X)  in  Bullrich's  notation  in  Eq.  (15)  corresponds  to  the  aerosol  phase  function  PA(50°). 
Now,  PA  =  xpR,  where  x  is  the  turbidity  determined  from  the  50°  scattering  intensity  and  pR  is  the 
Rayleigh  scattering  function.  (The  factor  2  was  dropped  because  the  measured  intensity  is  the  average 
of  the  polarized  components).  Assuming  that  the  aerosol  and  Ravleigh  phase  functions  are 
approximately  equal  near  50°  scattering  angle  for  clear  air,  PA(50J)  is  simply  a  relative  factor  that 
normalizes  the  aerosol  phase  function  to  the  Rayleigh  value  of  1 .06  at  50°,  and  renders  the  number  of 
particles  Independent  of  scattering  angle. 

The  exponent  v,  at  a  particular  altitude,  is  determined  from  the  wavelength  dependence  of  pA  (via 
x)  using  the  Angstrom  formula,  mentioned  previously.  That  is, 

PA  =  constant  At®  and  PR  ~  A.-4 

so  that  x  -  X,-a  +  41  =  X,e  '  v)  since  a  =  v  -  2  and  2  <  v  <  6.  It  follows  then  that  the  turbidity  has  a  positive 
wavelength  exponent,  a. 

x  ~  X“  where  a'  =  6  -  v  and  therefore  0  <  a’  <  4 
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The  exponent  v  In  the  particle  size  distribution  can  thus  be  estimated  from  the  quotient  of  turbidities 
measured  at  X[  and  Xj.-  That  is, 

li  v,  R  /]ogJxl/t2)\ 

t2  "  (xj  v“6_(log(M/A2)J  • 

The  angular  scattering  phase  functions  for  commonly  used  exponential  size  distribution  have  been 
computed  from  Mle  Theory13.  Consequently,  aerosol  particle  concentrations  that  give  rise  to 
measured  values  of  (JA  can  be  estimated  by  using  Eq.  (39). 

Such  a  procedure  is  deemed  reasonable  because  as  Harris  and  McCormick23  have  demonstrated, 
for  a  specific  range  of  particle  radii,  a  power-law  size  distribution  can  be  approximated  by  a  log¬ 
normal  distribution;  (shown  In  Section  3.2).  Then  too,  an  exponential  wavelength  dependence 
reflected  in  the  scattering  intensity  implies  that  a  power-law  size  distribution  can  be  assumed  without 
loss  of  generality  because  the  angular  scattering  factor,  although  a  function  of  size  parameter,  is  not 
an  explicit  function  of  wavelength.  The  latter  assertion  was  used  as  an  assumption  previously,  but 
Eq.  (39)  clearly  demonstrates  from  a  dimensional  perspective  that  PA(<f>)  is  independent  of  X.  As  a  first 
approximation  then,  the  aerosol  particle  concentration  at  representative  altitudes  can  be  inferred 
from  the  turbidity  data.  In  this  connection,  the  50°  normalization  factors.  PA(50°),  were  determined  by 
extrapolation  from  aerosol  phase  functions  tabulated  by  Bullrich  for  a  power-law  size  distribution  of 
radii  in  the  range  (0.04  pm  <  r  <  10  pm).  These  were  deemed  suitable  for  order-of-magnltude  estimates 
and  are  listed  in  Table  2.  along  with  the  foreward-to-backscatter  ratio.  However,  the  latter  ratio  is 
dependent  upon  the  upper  and  lower  limits  of  radii  used  in  the  model  and  therefore  varies  with  this 
range  interval,  which  is  an  unknown  in  the  nephelometer  data.  Specifically,  a  given  exponent  v  does 
not  uniquely  define  the  phase  function  and  a  strict  comparison  between  the  experimentally  observed 
ratios  and  those  in  Table  2  cannot  be  made. 


Table  2.  Theoretical  Scattering  Parameters  vs  Power  Law  Exponent 


V 

PA(50°) 

PA(39°)/PA(  1 50°) 

2.5 

2.3 

13.95 

3.0 

1.2 

13.88 

3.5 

0.78 

12.52 

4.0 

0.52 

10.04 

Figure  7a  is  an  altitude  profile  of  number  of  particles  per  cm3  greater  than  0. 15  pm  in  radius,  resulting 
from  the  conversion  of  turbidity  at  0.660  pm  to  number  density  by  using  Eq.  (39).  Table  1  lists  these 
data.  To  illustrate  a  typical  calculation,  the  turbidity  ratio  (for  0.660  pm  and  0.475  pm)  at  9  km  is 
about  3,  which  in  turn  implies  PA-X'0-6  orv  =  2.6. 

Additionally,  the  observed  value  of  x  (at  0.660  pm)  is  0.3  and  PR  =  1.84  x  10-9  cm-1  sr_1.  When 
these  values  are  used  in  Eq.  (39)  we  obtain  at  9  km, 


23.  Harris,  F.S.  and  McCormick,  M.P.  (1972)  Mie  scattering  by  three  polydispersions,  J.  ColL  and 
Inter.  Set  29:536-545. 
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a.  NO.  PARTICLES/cm3 
WITH  r>0.15nm 


b.  AEROSOL  MIXING  RATIO 


1  2  3  0  2  4  6  8  10  12  14  16 

N  (  >  0.15  Jim)  cm'3  MIXING  RATIO  (nSpr) 


Figure  7.  (a)  Number  of  Particles  vs  Altitude,  (b)  Particle  Mixing 
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N  sr  0.3  — — rg —  (>  0.1 5  pm  radius) . 


Extensive  measurements  by  Hofmann  and  Rosen24  using  balloon-borne  particle  counters  provide 
perhaps  the  best  set  of  data  available  for  comparison.  In  particular,  their  findings  address  the 
question  of  the  aerosol  background  during  non-volcanic  periods  and  the  altitudinal  variability  in  the 
concentration  of  particles  greater  than  0. 15  pm  in  radius.  Their  data  are  presented  in  terms  of  the 
aerosol  mixing  ratio,  which  is  the  quotient  of  the  number  of  particles  to  air  mass-density,  and  is 
defined  by 


Aerosol  Mixing  Ratio 


particles/cm3  particles 

milligrams  of  air/cm3  ~  milligrams  ' 


i 


The  air  mass-density  is  given  in  the  tables  of  the  Standard  Atmosphere  1976.  The  relevance  of  the 
profile  of  the  aerosol  mixing  ratio  lies  in  its  comparison  with  the  the  optical  turbidity  which  is 
proportional  to  the  number  of  aerosol  particles  per  air  molecule.  When  the  number  concentrations 
that  were  inferred  from  the  turbidity  measurements  were  converted  to  aerosol  mixing  ratio.  Figure 
(7b)  was  obtained.  This  profile  compares  favorably  with  the  data  of  Hofmann  and  Rosen  obtained  in 
1979  over  Laramie,  Wyoming,  and  also  shown  in  Figure  (7b).  Specifically,  they  found  peak  mixing 
ratios  between  10  and  20  particles/milligrams  in  the  stratospheric  Junge  layer  between  18  and  23  km 
altitude  along  with  an  increase  in  relative  concentration  of  small  particles  above  25  km.  However, 
wide  variability  in  size  and  mixing  ratio  was  observed  in  the  troposphere  with  significant  seasonal 
dependence.  For  example,  in  the  spring  of  1978  when  the  peak  mixing  ratio  for  the  stratosphere  was 
about  8  particles/milligram,  the  ratio  10  km  varied  by  at  least  an  order  of  magnitude,  from  2  to  20 
particles/milligram.  This  contrast  with  the  autumn  values  of  about  0.5  to  2  in  the  9  to  10  km  level  and 
the  stratospheric  values  remained  virtually  unchanged.  And  again,  the  1978-79  period  was  one  in 
which  volcanic  dust  influences  were  noticeably  absent  from  their  data.  The  high  degree  of  variability 
in  the  tropospheric  aerosols  appears  to  be  consistent  with  the  nephelometer  data  presented  in  this 
report. 

As  noted  above,  the  particle  concentrations  displayed  in  Figure  (7a)  refer  to  particles  greater 
than  0. 15  pm  for  an  assumed  power-law  exponent  in  the  size  distribution.  The  dissymmetry 
parameter,  on  the  other  hand,  is  a  measure  of  the  effect  of  all  particles  encountered  and  without 
additional  information  related  to  particle  size,  could  not  be  expected  to  provide  strict  comparison 
with  the  theoretical  model  used  in  determining  particle  number  densities.  From  a  qualitative 
viewpoint,  however,  the  compatibility  of  these  independent  sets  of  data  in  presenting  signatures  for 
the  presence  of  aerosols  is  evident.  Within  the  margins  of  experimental  error,  the  angular  scattering 
as  depicted  in  the  R-parameter  gives  a  reasonable  picture  of  the  change  in  the  aerosol  phase  function 
with  altitude  in  that  the  ratios  are  less  limited  to  values  characteristic  of  particles  with  radii  less  than 
about  0.2  pm.  This  can  be  ascertained  from  Van  de  Hulst10  wherein  the  forward-to-backward 
scattering  by  a  sphere  is  tabulated  as  a  function  of  size  parameter  and  the  magnitude  of  the  forward 
enhancement  as  a  function  of  particle  size  can  be  estimated.  His  Table  5  Section  7.22,  indicates  that 


24.  Hofmann,  D.F.  and  Rosen,  J.M.  (1981)  On  the  background  stratospheric  aerosol  layer,  J.  Atmos. 
Set  38:168-181. 
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for  size  parameters  less  than  2,  the  total  scattering  In  the  forward  hemisphere  exceeds  that  in  the 
backward  hemisphere  by  a  factor  less  than  10.70,  but  rapidly  increases  to  a  factor  of  62.0  when  the  size 
parameter  reaches  3.  Similarly,  a  size  distribution  such  as  Deirmendjian's  continental  haze  model25 
results  in  a  factor  16.5  as  the  ratio  between  the  value  of  phase  functions  at  30°  and  1 50°.  Other  models 
of  atmospheric  aerosols  such  as  Eltcrman22,  for  visibility  conditions  of  at  least  30  km,  also  indicate 
this  factor  lies  in  a  range  from  10-20,  depending  upon  the  relative  weighting  of  particles  with  radii 
greater  than  0.2  pm.  The  inferred  phase  functions  in  Table  1  are  reasonably  close  to  these  models. 
Notably,  the  relative  change  in  PA(150°)  is  more  pronounced  at  the  stratospheric  altitude  of  18  km 
than  in  the  troposphere,  which  is  consistent  with  an  increase  in  particle  size  in  the  Junge  layer.  The 
LOWTRAN  aerosol  phase  functions26  illustrate  this  behavior  of  the  backscatter,  for  a  wide  variety  of 
atmospheric  aerosol  models.  Therefore,  the  absolute  precision  in  values  for  the  dissymmetry  does  not 
preclude  an  informative  correlation  of  this  parameter  with  substantive  angular  scattering  properties 
of  aerosols,  and  reasonable  inferences  relating  to  particle  size  can  be  made  from  such  data. 

The  dissymmetry  parameters  for  the  1970  data  are  displayed  in  Figure  8  from  ground-level  to  15  km 
altitude  (due  to  failure  of  the  backscatter  photometer  at  15  km).  Those  results  are  consistent  with  the 
1976  data  in  which  a  comparable  wavelength  relationship  existed,  namely  dispersion  toward  the 
longer  wavelengths.  However,  in  Figure  8,  the  R-values  for  the  0.475  pm  wavelength  remained  less 
than  about  1.45  in  the  troposphere  (from  5-15  km)  indicating  virtually  constant  symmetry;  a  property 
of  Rayleigh  scattering.  Such  symmetry  was  absent  at  0.660  pm,  thereby  strengthening  the  suggestion 
made  in  Section  3.3  that  this  wavelength  was  insensitive  to  the  small  particle  component  in  the 
aerosol  distribution.  On  the  contrary,  the  1976  data  in  Figure  6  do  not  show  such  constant  R  values 
for  the  blue  wavelength.  Apparently,  both  the  subtle  differences  and  distinct  similarities  between  the 
sets  of  data  reveal  physically  significant  aspects  of  angular  scattering  from  atmospheric  particles 
encountered.  Specifically,  the  disparity  in  the  magnitude  of  variability  in  the  vertical  profiles 
indicates  the  contrasting  size  ranges  of  tropospheric  particles,  despite  the  compatibility  in 
wavelength  displacement. 

In  the  next  section,  a  connection  between  the  dissymmetry  and  polarization  is  made,  which 
correlates  the  primary  scattering  characteristics  of  small  particles.  This  correlation  provides  a  more 
quantitative  resolution  of  the  problem  of  the  "unusual"  turbidity/wavelength  behavior  evident  in  the 
1970  data.  In  particular,  the  polarization  of  the  scattered  light  can  be  used  to  separate  the  small 
particle  component  from  the  total  scattering  Intensity.  This,  in  turn,  facilitates  a  comparison  of  the 
size  distribution  of  large  particles  observed  for  the  two  contrasting  sets  of  data. 


3.5  Inferences  from  the  Polarizing  Properties  of  Small  Particles 

Perhaps  the  most  interesting  data  from  the  1970  balloon  flight  was  the  polarization  of  the 
scattered  light  at  the  90°  scattering  angle.  In  addition  to  angular  symmetry,  light  scattered  by  small 
particles  is  polarized  to  a  greater  degree  than  that  by  large  particles.  Therefore,  the  observed 


25.  Deirmendjian,  D.  (1969)  Electromagnetic  Scattering  in  Polydispersions,  Elsevier,  New  York. 

26.  Kneizys,  F.X.  et.  al.  (1983)  Atmospheric  Transmittance /Radiance:  Computer  Code  LOWTRAN  6, 

AFGL-TR-83-0187,  ADA137786. 
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ALTITUDE  (km) 


FORWARD-TO-BACKSCATTER  RATIO 


Figure  8.  The  Ratio  of  Forward-to-Backward  Scattering  (Dissymmetry)  vs  Altitude  (November  1970) 
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polarization  ratio  at  0.475  pm.  Figure  5,  can  be  used  to  examine  this  behavior  in  the  troposphere. 
Call27  suggested  a  procedure  for  investigating  the  relationship  between  polarization  and  aerosol 
particle  size,  using  a  LIDAR  system.  His  technique  was  adapted  to  the  nephelometer  observations  in 
the  following  manner:  For  a  particular  photometer,  the  photomultiplier  (PM)  detector  response  is 
linearly  related  to  the  volume  scattering  function  and  can  be  expressed  as, 

W(<)>s)  =  K  (4>s)  (41) 

where  K  is  the  calibration  factor  for  a  photometer  at  scattering  angle  <j>8  with  spectral  filter  of 
wavelength,  X.  From  Eq.  (2)  this  can  be  written 

W(4>s)  =  KloRX  PR(<f8)+oAxPA^s)  +  ^rxPR(Os)l 

(42) 


where  opX  is  the  scattering  coefficient  due  to  small  particles,  which  are  also  assumed  to  have  a 
Rayleigh  scattering  phase  function  PR(<j>s).  The  other  terms  were  defined  previously  in  connection  with 
Eq.  (2).  The  PM  tube  is  insensitive  to  the  polarization  of  the  scattered  light,  but  the  vertical  and 
horizontal  polarized  intensities  can  be  detected  sequentially  with  properly  oriented  polarizers.  The 
balloon-borne  photometers  used  this  arrangement.  The  detector  responses  then  take  the  form 

=  CWV(<J>3)  =  K1(aR^PRV(<jts)  +  ^axPav^s)  +  ^pxPrv^s^  ( 43) 

and 

IH(«t»J  =  CWH(*S)  =  Kjl0R^PRH(4»s)  +  OaxPah^s)  +  CTpi  PkhI<I>sJ1  (44) 

where  K,  is  a  constant  which  includes  the  linear  relationship  between  the  PM  tube  response  and  the 
Intensity.  The  subscripts  V  and  H  on  the  respective  phase  functions  simply  refer  to  the  vertical  and 
horizontal  components.  Note  that  the  scattering  coefficients  o  are  independent  of  scattering  angle. 

The  ratio  of  Eq.  (43)  to  Eq.  (44)  gives  a  quantity  with  the  calibration  factor  eliminated,  and  is  simply 
the  polarization  ratio, 

ly(^3)  _  CTRi.  Prv^s)  +  gAXPAy(^8^  +  °pX  Prv^s)  .  . 

Ih<*s)  "  <*Rt  PRH(<t>s)  +  <*axPah(4>s)  +  °PX  Prh(<I>s) 

At  <)>s  =  90°  this  expression  is  simplified,  because  PRH  (90°)  is  proportional  to  cos2  <ps,  which  equals  zero 
at  <t>„  =  90°.  Therefore  for  the  photometer  at  <)>s  =  90°  we  obtain 

(y  _  ctrhPrv(90°) 

“  Wah^) 


Pav(90°) 
d90°) 


■ah' 


gpX  PRy(90°) 

^axPahOO0) 


(46) 


Using  the  relationship 

qPx  _  <*pX  gRX 
ctrX  ctax 


27.  Call,  R  (1967)  Measurement  of  Atmospheric  Aerosols  by  Polarized-Laser  Ught  Scattering,  Unlv. 
Microfilms,  Inc.,  Ann  Arbor,  MI. 
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we  get 


ly  _  °rX  f .  CTPXA  Prv(90°)  PAy(90°) 

I«  '  <*Ax  +  <TrxJpah<90°)  +  PAHO0°)  * 


Finally, 


(1  +  opX/oRX) 


Iy/IH-  PAv(90°)/Pah(90o)  oaX 
Prv190°) /Pah(90°)  orX  • 


(48) 


(49) 


This  result  shows  that  the  small  particle  turbidity  factor,  (1  +  ap/aR)  can  be  determined  from  the 
polarization  ratio  and  is  directly  related  to  the  turbidity,  oA/aR.  It  is  important  to  note  that  at 
X  =  0.475  pm,  this  value  of  aA/oR  arises  from  both  small  particles  and  large  particles  and  can  be 
equated  to  the  observed  value  at  the  45°  scattering  angle,  since  a  is  independent  of  scattering  angles. 
For  example  if  oP  =  0;  that  is,  there  are  no  small  particles,  the  turbidity  due  to  large  particles  alone  is 
from  Eq.  (49) 


Prv/P 


/_  /_  ,  _  -  RV'  AH 

1°a/CTr)^p=0  "  (Iv/Ih)o  -  PAv/P/ 


AH 


(50) 


The  observed  turbidity  is  therefore 

.  ,  ,  ,  ,  v  M  ,  ,  (Iv/Ih^O  ~  (PAy/PAH^ 

(CTA/CTR)0bs  “  (CTA/aR>np=o  (1  +  °p/aR)  (iv/iH)  _  (PAV/PAH)  • 

It  is  thus  evident  that  the  observed  turbidity  is  greater  than  that  due  to  large  particles  alone  when 
( 1  j  X  ly/lw  ~  PAy/PAH 

(1  +aP/oR)  >  (Iv/Ih,0_  Pav/pah 


(51) 


(52) 


or 


<rP/cR  > 


(Iv/Ih)o~  PAv/PAh 


(53) 


A  further  implication  from  Eq.  (50)  and  this  inequality  is  that 


(°a/ctr)(jp=o  < 


(gP/qR)  Prv/PAh 


(54) 


It  is  evident  therefore,  that  the  aerosol  turbidity  factor  (aAX/oRX)  in  Eq.  (49)  is  a  function  of  the  relative 
concentrations  of  small  particles  which  cause  the  enhanced  polarization  ratio. 

Inspection  of  Figure  5  indicates  that  near  ground  level  and  in  the  stratosphere  (altitude  >16  km), 
both  regions  where  large  particles  were  dominant,  Iv/Ih  18  (ess  than  6.  On  the  other  hand,  this  ratio 
reaches  values  between  8  and  10  in  the  troposphere,  where  the  small  particle  contribution  to  both  the 
polarization  and  turbidity  must  be  considered.  Now,  Eq.  (49)  provides  a  means  to  determine  the  factor 
(1  +  oP/oR),  at  0.475  pm,  from  the  observables  lv/ln  and  (oA  /oR)  if  suitable  aerosol  phase  functions  can 
be  found.  Particularly  applicable  in  this  connection  were  phase  functions  used  by  Call27  in  his  LIDAR 
work,  because  he  derived  and  plotted  the  two  polarized  components  as  a  function  of  scattering  angle. 
He  assumed  an  exponential  size  distribution  and  refractive  index  of  1.5  for  the  aerosol  particles.  The 
pertinent  parameters,  giving  the  best  fit  to  the  measured  data  are  given  in  Table  3. 
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Table  3.  Theoretical  Polarization  Parameter  vs  Power-Law  Exponent 


V 

Size  Range 

Pr^J/P^OO0) 

PAv(90°)/Pah(90°) 

3.5 

.1  <  r<  1  pm 

7.93 

2.35 

3.0 

.1  <r<  1  pm 

7.53 

1.84 

2.5 

.1  <r<  1  pm 

7.00 

1.35 

To  identify  the  proper  phase  function  at  a  particular  altitude,  the  dissymmetry  parameter 
provides  a  source  for  testing  the  data  consistency  because  it  also  contains  the  small  particle 
parameter,  (1  +  oP/oR).  This  condition  arises  in  a  manner  similar  to  Eq.  (6)  when  the  small  particles 
are  included.  That  is,  in  view  of  the  assumption  that  the  small  particles  are  Rayleigh  scatterers  and 
thus  have  the  same  phase  function  as  air  molecules,  the  ratio  of  forward-to-backward  scattering 
intensities  is  now 

R  _  (1  +  qP/gR)  Pr(45°)  +  (qA/qR)  PA(45°) 

'  (1  +  oP/oR)  Pr(135°)  +  (oA/oR)  Pa(135°)  U 


where  the  symmetric  scattering  angles  for  the  1970  data  were  45°  and  135°.  In  this  case  the  ( 1  +  cos2  <j>s) 
Rayleigh  factor  gives  Pr(45°)  =  Pr(135°)  =1.12  and  Eq.  (55)  becomes 

1.12  (1  +  qp/aR)  +  tPa(45°) 

1. 12  (1  +  oP/qR)  +  xPA(135°)  1 


and  x  =  (qA/qR )  is  the  turbidity  due  to  large  (Mie)  particles  alone  as  observed  with  the  photometer  at 
45°.  Rearranging  this  expression  gives 

Pa(45°)  =  1.12  (1  +  Op/OR)  f~ •)  +RPa(135°)  (57) 


and  with  the  assumption  that  the  phase  function  is  independent  of  wavelength,  used  previously  in 
Section  3.4,  we  arrive  at  an  expression  analogous  to  Eq.  (37) 

+  1581 

This  relation  assumes  that  X;  is  due  to  only  Mie  particles.  In  the  interpretation  of  the  1970  turbidity 
data,  (Figure  3),  it  was  asserted  that  the  observed  values  of  x2  at  0.475  pm  were  masked  by  the  small 
particles  in  the  troposphere.  To  apply  Eq.  (58)  it  is  necessary  to  infer  the  probable  values  of  x2  from 
those  at  X.  =  0.660  pm,  by  assuming  a  turbidity-wavelength  dependence  consistent  with  the 
polarization  data.  That  is,  the  choice  of  aerosol  phase  function  that  best  matches  the  polarization  and 
turbidity  parameters  in  Eq.  (49)  establishes  the  exponent  v  in  the  size  distribution  as  well  as  the  X- 
dependence  in  the  turbidity  from  the  relation 

x-X6 ‘v  and  thus  ^  .  (59) 

Equation  (58)  can  finally  be  expressed  in  terms  of  the  v  and  Xj,  the  observed  turbidity  at  X  =  0.660  pm. 
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(60) 


pa'135°>  -  EFTS  11  %y°pl  [».  -  « (^-J  '  "  - • 

Again,  this  expression  Is  quite  sensitive  to  the  exponent  v,  which  cannot  be  chosen  arbitrarily  if  the 
data  are  to  be  consistent  and  because 


Xt  \6-v  R^l 

h)  > 


(61) 


If  the  Pa(135°)  has  positive  values. 

This  procedure  was  employed  to  examine  the  variability  In  aerosol  phase  functions  up  to  only  15 
km  for  the  1970  data  because  of  the  failure  of  the  135°  photometer  previously  mentioned. 

Moreover,  Initial  guidance  for  making  judicious  selection  of  values  of  v  comes  from  the  fact  that  the 
values  of  x  and  R  at  0.660  pm  were  comparable  for  both  flights  in  1970  and  1976.  This  suggested  that 
the  size  distributions  of  large  particles  were  not  appreciably  different  In  the  troposphere.  The 
analytical  procedures  thus  Involved  three  steps.  First,  an  estimate  of  the  v-value  satisfying  condition 
(61)  was  made.  Second,  this  value  defined  the  phase  function  used  In  Eq.  (49)  at  a  particular  altitude,  to 
determine  (1  +  oP/aR)  from  the  observables  IV/IH  and  (cA/oR).  Third,  the  value  so  determined  was  used  to 
find  PA(135°)  from  Eq.  (60)  and  subsequently  PA(45°)  from  Eq.  (57).  Finally,  In  the  same  manner  as  in 
Section  3.4,  it  was  possible  to  estimate  the  number  of  particles  greater  than  0.15  pm  for  a  given  v-value 
from  the  observed  turbidity  values  at  0.660  pm.  Table  4  lists  the  values  of  the  above  parameters  for  the 
1970  flight.  Included  In  the  table  Is  the  small  particles  turbidity.  xP  =  oP/aR  and  two  other  quantities; 
the  number  of  small  particles  (radius  >  0.04  pm)  per  cubic  centimeter,  and  the  small  particle  mixing 
ratio.  The  number  concentration  of  small  particles  can  be  estimated  from  Eq.  (39)  by  substituting  the 
parameter  values;  Pp  =  xP  Pn.  X  =  0.475  pm,  r  =  0.04  pm,  v  =  6  and  PR  (45°)  =1.12.  As  for 


Table  4.  Scattering  Properties  and  Particle  Number  Concentrations  (November  1970) 


Altitude 

(km) 

R1 

D 

Ra 

a 

KQj 

Pa(45”) 

Pa(135") 

a 

N>0.15  mm 

N>0.04  mm 

1  - 

on2 

mg 

mg 

2.20 

1.9 

1.50 

0.71 

1.10 

2.04 

0.25 

8.16 

2.68 

3.0 

3.50 

3.98 

1800 

2050 

2.15 

1.8 

1.45 

0.63 

1.22 

2.19 

0.28 

7.82 

2.87 

2.8 

3.20 

4.03 

1837 

2312 

1.70 

1.0 

1.25 

0.40 

1.40 

2.46 

0.44 

5.59 

3.17 

2.5 

2.18 

3.04 

1792 

2503 

2.15 

1.2 

1.35 

0.48 

1.24 

3.12 

0.33 

9.45 

3.17 

2.5 

1.40 

2.31 

1426 

2214 

1.70 

1.1 

1.25 

0.44 

1.40 

2.46 

0.44 

5.59 

3.17 

2.5 

1.59 

2.74 

1477 

2542 

8 

1.68 

0.9 

1.25 

0.35 

0.80 

2.11 

0.40 

5.27 

3.06 

2.6 

1.23 

2.35 

740 

1378 

9 

1.60 

0.7 

1.20 

0.28 

0.85 

2.18 

0.25 

8.72 

3.17 

2.5 

0.87 

1.86 

715 

1481 

10 

1.50 

0.8 

1.20 

0.30 

0.63 

1.49 

0.23 

6.48 

2.65 

3.0 

0.69 

1.64 

459 

1093 

11 

1.76 

1.2 

1.30 

0.39 

0.80 

2.30 

0.58 

3.96 

3.06 

2.6 

1.04 

2.78 

513 

1372 

12 

1.70 

1.3 

1.35 

0.49 

0.80 

2.29 

0.71 

3.22 

2.65 

2.8 

0.84 

2.53 

442 

1331 

13 

1.50 

1.2 

1.20 

0.38 

0.74 

1.91 

0.73 

2.61 

3.17 

2.5 

1.17 

3.98 

344 

1170 

14 

1.75 

1.2 

1.40 

0.53 

0.46 

2.14 

0.63 

3.40 

2.28 

3.5 

0.41 

1.58 

187 

122 

15 

1.95 

1.3 

1.40 

0.49 

0.19 

1.46 

0.20 

6.80 

2.65 

3.0 

0.78 

3.43 

60 
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large  particles,  the  small  particle  mixing  ratio  is  simply  the  quotient  of  the  number  concentration 
divided  by  the  air  mass  density. 
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; 


) 


It  Is  readily  apparent  from  Table  4  that  the  number  concentration  of  small  particles  exceeds  the 
large  particle  concentration  by  a  factor  of  at  least  0.5  x  103,  In  the  altitude  range  5-12  km.  This  result 
Is  consistent  with  the  qualitative  estimates  deduced  from  Eq.  (25)  in  Section  3.3,  when  allowance  is 
made  for  the  v-values,  indicated  for  the  large  particles.  That  is,  the  latter  values  are  less  than  3,  while 
Eq.  (25)  assumed  v  =  4. 

In  the  8-15  km  region,  the  large  particle  concentrations  for  the  two  flights  are  closely 
comparable.  (No  data  were  used  below  8  km  for  the  1976  flight,  because  the  balloon  ascent  rate  was  too 
rapid.) 

This  region  Is  particularly  significant,  however,  because  It  has  been  previously  Identified  as  one 
in  which  backscaiter  from  ground  based  LIDAR  systems  reached  minimum  values,  and  was  thus 
considered  appropriate  for  data  normalization  purposes28. 

4.  CONCLUSIONS 

The  nephelometer  scattering  observations  provide  self-consistent  vertical  profiles  of  aerosol 
scattering  effects  In  the  troposphere  and  stratosphere.  The  optical  properties  of  the  stratosphere  over 
New  Mexico  showed  little  change  from  1970  to  1976  even  though  the  Intervening  period  included  a 
major  volcanic  dust  event  and  the  rare  observation  of  a  nacreous  cloud  at  25  km6.  By  contrast, 
however,  the  troposphere  exhibited  wide  variability.  Including  evidence  of  optical  scattering  effects 
attributable  to  the  gamut  of  aerosol  particles  from  condensation  nuclei  to  Mie  particles  even  in 
presumably  clear  air.  In  particular  the  data  suggest  that  small  particles  that  exhibit  Rayleigh 
scattering  cannot  be  neglected  In  a  realistic  picture  of  aerosol  optical  effects  in  the  atmosphere.  That 
such  particles  may  exist  in  large  concentrations  In  a  very  clear  sky  is  amply  demonstrated  by  the 
formation  of  jet  contrails,  and  Indeed  they  are  critical  to  the  growth  of  clouds,  as  nucleating  agents. 
The  discrepancy  between  the  results  obtained  with  LIDAR  probing  and  those  from  other  atmospheric 
optical  techniques,  as  pointed  out  by  Rozenberg29,  finds  further  confirmation  in  the  data  from  the 
nephelometer  flights.  Specifically,  the  values  of  optical  turbidity  observed  were  at  least  a  factor  of  2 
greater  than  those  determined  by  Ildar  during  non-volcanic  periods.  Additionally,  the  nephelometer 
data  showed  an  increase  in  turbidity  and  aerosol  mixing  ratio  above  23  km,  despite  a  reduction  in 
concentration  of  particles  with  radii  greater  than  0. 15  pm.  This  result  is  apparently  consistent  with 
the  finding  of  Hofmann  and  Rosen24  relating  to  increased  numbers  of  small  particles  at  altitudes 
above  the  Junge  layer;  a  result  not  evident  from  atmospheric  sounding  with  lidar,  but  clearly 
manifested  in  the  occurrence  of  high  altitude  clouds. 

The  nephelometer  observations  support  the  finding  by  others  indicating  the  existence  of  an  ever¬ 
present  background  of  aerosol  particles  in  the  earth’s  atmosphere  against  which  the  recent  dust 
Incursions  from  eruptions  such  as  Mt  St  Helens  and  El  Chichon  should  be  measured.  Morever,  this 
background  of  aerosols  has  been  shown  to  be  somewhat  self-sustaining  and  perhaps  natural  in  origin. 


28.  Russell,  P.B.  et.  al.  (1979)  Methodology  for  error  analysis  and  simulation  of  Lidar  aersosol 

measurements,  Appl  Opt 18:3783-3797. 

29.  Rozenberg,  G.V.  and  Sandormirsky,  A.B.  (1971)  The  optical  stratification  of  atmospheric 

aerosols,  fzu.  Atmospheric  and  Oceanic  Physics  7:737-749. 
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although  uncertainties  remain  as  to  a  clear  source,  because  of  anthropogenic  as  well  as  continuing 
volcanic  Influences. 


32 


References 


1.  Northam,  G.B.  et.  al.  (1974)  Dustsonde  and  Ildar  measurements  of  stratospheric  aerosols:  a 

comparison.  Appl.  Opt  13:2416-2421. 

2.  Russell,  P.B.  et.  al.  (1975)  Results  of  stratospheric  lidar  observations.  Fourth  Conference  on  the 

Climatic  Impact  Assessment  Program,  Cambridge.  MA,  4-7  February. 

3.  McCormick,  M.P.  and  Fuller.  W.H.  (1975)  Lidar  measurements  of  intense  stratospheric  dust 

layers,  Appl.  Opt.  14:4-5. 

4.  Barteneva,  O.D.  (1960)  Scattering  functions  of  light  in  the  atmospheric  boundary  layer.  Bull 

Acad.  Set  USSR  Geophys.  ScL,  1237-1244. 

5.  Gibson,  F.W.  and  Dearborn.  F.K.  (1971)  Atmospheric  Optics  Measurements  with  a  Balloon  Borne 

Nephelometer,  AFCRL-TR-7 1-0455.  AD736408. 

6.  Gibson,  F.W.  and  Volz,  F.E.  (1972)  High  altitude  measurements  of  the  optical-scattering 

properties  of  the  atmosphere,  presented  at  Optical  Society  of  America  Meeting.  New  York,  NY, 
April. 

7.  Gibson,  F.W.  ( 1 976)  In-situ  photometric  observations  of  angular  scattering  from  atmospheric 

aerosols,  Appl  Opt  15:2520-2533. 

8.  Meinel.  A.B.  and  Metnel,  M.P.  (1975)  A  stratospheric  dust  event  of  November  1974,  Science 

188:477-461. 

9.  Junge.  C.E.,  Chagnon,  C.W.,  and  Manson,  J.E.  (1961)  Stratospheric  aerosols.  J.  Meteorol  18: 

81-98. 

10.  Van  de  Hulst.  H.C.  ( 1957)  Light  Scattering  by  Small  Particles,  John  Wiley  &  Sons,  New  York. 

11.  Sassen,  K.  (1981)  Infrared  (10.6  pm)  scattering  and  extinction  in  laboratory  water  and  ice  clouds. 

Appl  Opt  20:185-193. 

12.  Kerker,  M.  (1969)  The  Scattering  of  Light  and  Other  Electromagnetic  Radiation,  Academic  Press. 

New  York. 

13.  Bullrich,  K.  (1964)  Scattered  radiation  in  the  atmosphere  and  the  natural  aerosol.  Advances  in 

Geophysics  10:101-257. 

14  Junge.  C.E.  (1961)  Vertical  profiles  of  condensation  nuclei  in  the  stratosphere,  J.  Meteorol. 
18:501-509. 


33 


15.  Went,  F.W.  (1964)  The  nature  of  Aitkin  condensation  nuclei  In  the  atmosphere.  Nature  51: 1259- 

1267. 

16.  Whitby,  K.T.  (1975)  Modeling  of  Atmospheric  Aerosol  Particle  Size  Distribution,  Prog.  Rep.  253, 

Particle  Technology  Laboratory,  Univ.  of  Minnesota,  Minneapolis. 

17.  Junge,  C.F.  and  McLaren,  E.  (1971)  Relationship  of  cloud  nuclei  spectra  to  aerosol  size 

distribution  and  composition,  J.  Atmos.  Set  28:382-390. 

18.  Nilsson,  B.  (1979)  Meteorological  influence  on  aerosol  extinction  in  the  0.2  pm  -  40  pm 

wavelength  range,  Appi  Opt  18:3457-3473. 

19.  Shettle,  E.P.  and  Fenn,  R  W.  (1979)  Models  for  the  Aerosol  of  the  Lower  Atmosphere  and  the 

Effects  of  Humidity  Variations  on  their  Optical  Properties,  AFGL-TR- 79-0214,  ADA085951. 

20.  Rossler,  F.  (1972)  Aerosol  layers  in  the  atmosphere,  Space  Research  XII  29:423-431. 

21.  Pilipowsky,  S.  et.  al.  (1968)  Investigation  of  the  stratospheric  aerosol  by  infrared  and  Lidar 

techniques,  J.  Geophys.  Res.  73:7553-7560. 

22.  Elterman,  L.  (1966)  Aerosol  measurements  in  the  troposphere  and  stratosphere.  Appl  Opt 

5:1769-1775. 

23.  Harris.  F.S.  and  McCormick,  M.P.  (1972)  Mie  scattering  by  three  polydispersions,  J.  Coll,  and 

Inter.  Set  29:536-545. 

24.  Hofmann.  D.F.  and  Rosen.  J.M.  (1981)  On  the  background  stratospheric  aerosol  layer.  J.  Atmos. 

Set  38:168-181. 

25.  Deirmendjian,  D.  (1969)  Electromagnetic  Scattering  in  Poly  dispersions.  Elsevier,  New  York. 

26.  Kneizys,  F.X.  et.  al.  (1983)  Atmospheric  Transmittance /Radiance:  Computer  Code  LOWTRAN  6. 

AFGL-TR-83-0187,  ADA137786. 

27.  Call.  R  (1967)  Measurement  of  Atmospheric  Aerosols  by  Polartzed-Laser  Light  Scattering.  Univ. 

Microfilms,  Inc.,  Ann  Arbor,  MI. 

28.  Russell,  P.B.  et.  al.  (1979)  Methodology  for  error  analysis  and  simulation  of  Lidar  aersosol 

measurements.  AppL  Opt. 18:3783-3797. 

29.  Rozenberg,  G.V.  and  Sandormirsky,  A.B.  (1971)  The  optical  stratification  of  atmospheric 

aerosols,  Izu.  Atmospheric  and  Oceanic  Physics  7:737-749. 


34 


